基于贝叶斯框架的各向异性页岩储层岩石物理反演技术   总被引：1，自引：0，他引：1       下载免费PDF全文

刘财  符伟  郭智奇  刘喜武  刘宇巍 《地球物理学报》2018,61(6):2589-2600

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究，针对页岩黏土含量高、层间微裂缝发育等特点，利用Backus平均理论描述页岩黏土矿物弹性参数，利用Chapman理论计算与水平微裂缝有关的VTI各向异性，并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法，并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息，并以测井纵、横波速度作为约束条件建立反演的目标函数，同时利用粒子群算法进行最优化搜索.计算结果表明，基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高，并且在一定范围内变化，这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关，同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关，而与横波各向异性参数γ相关性较小.另外，页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.  相似文献   

A model‐based data‐driven dictionary learning for seismic data representation          下载免费PDF全文

Can Evren Yarman  Rajiv Kumar  James Rickett 《Geophysical Prospecting》2018,66(1):98-123

Planar waves events recorded in a seismic array can be represented as lines in the Fourier domain. However, in the real world, seismic events usually have curvature or amplitude variability, which means that their Fourier transforms are no longer strictly linear but rather occupy conic regions of the Fourier domain that are narrow at low frequencies but broaden at high frequencies where the effect of curvature becomes more pronounced. One can consider these regions as localised “signal cones”. In this work, we consider a space–time variable signal cone to model the seismic data. The variability of the signal cone is obtained through scaling, slanting, and translation of the kernel for cone‐limited (C‐limited) functions (functions whose Fourier transform lives within a cone) or C‐Gaussian function (a multivariate function whose Fourier transform decays exponentially with respect to slowness and frequency), which constitutes our dictionary. We find a discrete number of scaling, slanting, and translation parameters from a continuum by optimally matching the data. This is a non‐linear optimisation problem, which we address by a fixed‐point method that utilises a variable projection method with ?₁ constraints on the linear parameters and bound constraints on the non‐linear parameters. We observe that slow decay and oscillatory behaviour of the kernel for C‐limited functions constitute bottlenecks for the optimisation problem, which we partially overcome by the C‐Gaussian function. We demonstrate our method through an interpolation example. We present the interpolation result using the estimated parameters obtained from the proposed method and compare it with those obtained using sparsity‐promoting curvelet decomposition, matching pursuit Fourier interpolation, and sparsity‐promoting plane‐wave decomposition methods.  相似文献   

Gully head modelling: A Mediterranean badland case study          下载免费PDF全文

Dino Torri  Jean Poesen  Mauro Rossi  Valerio Amici  Daniele Spennacchi  Cati Cremer 《地球表面变化过程与地形》2018,43(12):2547-2561

Predicting the location of gully heads in various environments is an important step towards predicting gully erosion rates. So far, field data collection and modelling of topographic thresholds for gully head development has mainly focused on gullies that formed in forested areas, rangelands, pastures and cropland. Such information for gullies in badlands however is very scarce. Therefore, this paper aims to extend the database on gully head topographical thresholds through data collection in a badland area and to improve the prediction of gully heads forming at sites with a very low erosion resistance value. For this, we chose a badland site located in central Italy that is characterized by biancana forms and both active and dormant gullies. The definition of the conditions under which present‐day gully heads developed allowed a better modelling of the gully head threshold equation, with modification of a previous model and the exemplification of how to use the updated model. The model shows that the resistance to gully head retreat depends on slope gradient and drainage area at gully heads, land use at the moment of gully development (as numerically expressed using parameters derived from the Runoff Curve Number method), surface rock fragment cover, presence of joints, pipes, and factors/processes affecting detachment rate. This study attempted to better understand environmental conditions that control the development of gully heads in badlands through a combination of field data collection of gully heads, an analysis of land use changes over 10 centuries, focusing on the period 1820–2005, and land use management through repeat photography and a critical examination of historical documents. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.  相似文献   

近距离上层隧道施工对下层隧道的扰动影响          下载免费PDF全文

唐强森  刘涛 《地震工程学报》2018,40(3):591-596,603

针对近距离上层隧道施工时给下层隧道稳定性造成扰动影响的现状,基于简支土体梁力学模型,研究下层隧道顶部在周期性扰动力作用下的系统振动特性,得到顶部垂直位移的动态响应数学表达式,为定量掌握最大位移峰值的大小与出现时刻提供理论依据,指出考虑振动因素时上下层隧道安全距离的确定方法。结合具体算例探索扰动力等因素对最大位移的作用规律,结果表明:扰动力振幅和隧道跨度分别增加1倍和2倍时,最大位移峰值相应增大1倍和23.9倍;隧道层间垂距和土体弹性模量分别增加2倍和3倍时,最大位移峰值相应减小88.2%和70%。当主要影响因素变化时,最大位移的波动大周期延长、波动震荡性加剧,特别是峰值处的震荡幅值增大,隧道顶部稳定性受到明显影响。  相似文献   

广州市轨道交通二十一号线工程路基抗震性测试研究   总被引：1，自引：0，他引：1       下载免费PDF全文

李同安 《地震工程学报》2018,40(4):759-764

传统路基动模量法大都进行单路基抗震性能分析,缺乏复合路基抗震性能分析,存在较高的偏差。对广州市轨道交通二十一号线工程路基抗震性进行测试,分析该地地质条件后,通过标准贯入实验法判断原始路基土质液化情况,使用瑞典条分法测试原始路基的抗震性能,判别出该路基液化土层滑动安全系数低、抗震稳定性差,应对其进行复合路基加固。通过圆弧滑动法对加固后的复合路基中碎石桩桩网结构路基以及CFG桩桩网结构路基进行抗震性测试。实验结果表明,所提方法抗干扰性能强且测试结果精度高。  相似文献   

2013年岷县—漳县M_S6.6地震前水氡浓度的临界慢化现象研究     

王熠熙  李赫  王博  杨朋涛  王俊  向阳  王喜龙  李悦 《地震》2018,38(1):128-138

基于临界慢化理论, 以2013年7月22日甘肃岷县—漳县M_S6.6地震前南北构造带及邻区水氡浓度观测资料为例, 逐一计算了表征临界慢化现象的自相关系数和方差。 研究结果表明, 2013年岷县—漳县M_S6.6地震前, 10个台点水氡浓度资料表现出较明显的临界慢化现象, 且出现慢化现象的观测点空间分布具有一定丛集性, 慢化持续时间呈现出由南向北迁移的特征。 综合分析认为, 临界慢化方法可以有效识别水氡资料蕴含的慢化信息, 这些早期异常信号对判定前兆异常所处阶段以及深入理解前兆资料变化的物理内涵具有重要的科学意义。  相似文献   

华南陆缘高热流区的壳幔温度结构与动力学背景          下载免费PDF全文

张健  王蓓羽  唐显春  董淼  艾依飞 《地球物理学报》2018,61(10):3917-3932

华南陆缘是我国重要的矿产、地热资源区.晚中生代以来，在太平洋板块西向俯冲，地幔热对流活动共同作用下，该区出现多期岩浆-热事件和大规模爆发式成矿作用.在前人研究基础上，本文利用地表热流观测资料、地震剪切波资料、重力位球谐系数，计算了壳-幔温度结构，分析了动力学背景.计算结果表明：华南陆缘东南沿海地带，地壳10 km以浅温度达200℃以上，居里点温度475℃，莫霍面平均温度550℃.地壳浅层较热，花岗岩中放射性元素衰变放热是地壳浅层地下水热活动的重要热源，但地壳总体温度不高，为"冷壳热幔"型热结构.地幔中，90 km深度，温度950~1250℃；120 km深度，温度1050~1400℃；150 km深度，温度1200~1450℃；220 km深度，温度1500~1700℃."热"岩石圈底界深度在110~150 km之间，西深东浅.岩石圈内，地幔应力场为挤压-伸展相间格局；岩石圈之下，地幔应力场为一个以南昌为中心、长轴NE-SW向的椭圆.分析认为，晚中生代以来，太平洋板块的西向俯冲，导致华南陆缘在区域性SE向地幔对流背景上叠加局域性不稳定热扰动，在175~85Ma期间，上地幔物质向上流动，形成不同的岩浆活动高峰期.同时，岩石圈地幔受俯冲洋壳流体的影响，含水量高，黏度小，在地幔流切向应力场作用下，岩石圈底界由西向东"波浪"状减薄.现今岩石圈之下仍具备地幔小尺度热对流温度条件，但除地表浅层外，地壳整体温度不高，岩石圈构造稳定.  相似文献   

3D P-WAVE VELOCITY STRUCTURE AT THE NORTHEASTERN MARGIN OF ORDOS BLOCK     

HAN Xiao-ming  Liu Fang  ZHANG Fan  CHEN Li-feng  LI Juan  LI Shuan-hu  YANG Hong-ying 《地震地质》2018,40(1):215-231

Using the 7 100 absolute first arrivals of P waves and 91 513 relative P arrival times of 726 events at the northeastern margin of the Ordos block since 2009, the 3D fine structure of P wave velocity within the depth of 15km in the crust was inverted by the double difference seismic tomography method. The results show that there exist obvious high-speed continuous bodies in the northwest of the study area, and their lateral areas increase gradually with depth, while the velocity of east and south is relatively low. The velocity inhomogeneity exists and differs at different depths. The lateral differences of velocity are related to seismicity and faults. The 5~15km depth profile shows that earthquakes tend to occur in the area with relatively high velocity or high speed transition zones, which to some extent reflects the fragility of regional crustal media and the strong differential movement of faults in vertical and horizontal directions where the crust body is easy to absorb and store strain energy and generate major earthquakes. A "Y"-shape low-velocity channel is present in the lower crust around Liangcheng, corresponding to the NW-trending Heilaoyao-Shahukou fault set, which may reveal the migration path of the Late Tertiary-Quaternary basalt eruption. The Helingeer M6.2 earthquake in 1976 was related to the formation of the locking section of the thermal welding in this area. The three-dimensional fine structure of P wave velocity presented in this paper provides intuitive seismological evidence for physical and chemical properties of crustal media and the deep tectonic environment of earthquake preparation.  相似文献   

Focal Mechanism and Focal Depth of the May 22, 2016 MS4.6 Earthquake in Chaoyang, Liaoning     

Zhao Xing  Zhang Fan  Han Xiaoming  Wang Shubo  Wei Jianmin 《中国地震研究》2018,32(1):100-112

An earthquake with M_S4.6 occurred at 17:08 p.m., May 22, 2016 in Chaoyang County, Liaoning Province. We used the P-wave first motion method, TDMT method, and CAP method to determine the focal mechanisms and the PTD method and sPn-Pn method to determine the focal depth. The focal mechanism results of the three methods are consistent. The depth results of the CAP method, PTD method and sPn-Pn method are close. We used the double difference location method to relocate earthquakes in 2009-2016, and obtained the strikes and dip angles of the small earthquake distributions with the help of simulated annealing algorithm and gauss Newton algorithm fitting. According to the focal mechanism results, the depth results, the characteristics of small earthquake distributions and the structural characteristics of the source area, the seismogenic fault strike is NEE and the main pressure force direction is NNW. The earthquake focal mechanism is for a normal fault type with a little left-lateral strike slip motion.  相似文献   

The NE Directed Seismicity Belt in Tibet after the MS8.1 Nepal Earthquake and Its Predictive Significance     

Wang Shuangxu  Zhu Liangyu  Xu Jing  Ji Lingyun  Jiang Fengyun 《中国地震研究》2018,32(1):119-129

After the 2015 M_S8.1 Nepal earthquake, a strong and moderate seismicity belt has formed in Tibet gradually spreading along the northeast direction. In this paper, we attempt to summarize the features and investigate the primary mechanism of this behavior of seismic activity, using a 2-D finite element numerical model with tectonic dynamic settings and GPS horizontal displacements as the constraints. In addition, compared with the NE-trending seismicity belt triggered by the 1996 Xiatongmoin earthquake, we discuss the future earthquake hazard in and around Tibet. Our results show that:the NE-directed seismicity belt is the response of enhanced loading on the anisotropic Qinghai-Tibetan plateau from the Indian plate and earthquake thrusting. Also, this possibly implies that a forthcoming strong earthquake may fill in the gaps in the NE-directed seismicity belt or enhance the seismic hazard in the eastern (the north-south seismic zone) and western (Tianshan tectonic region) parts near the NE-directed belt.  相似文献